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I. INTRODUCTION 

Nanoscale semiconductor heterostructures which are 
under intensive studies for more than 40 yearsi and par- 
ticularly quantum dots (QDs) play an important role 
in modern physics. Electrons and holes in such struc- 
tures are confined in all directions and their states be- 
come quantized. The interaction of the light with QDs 
can be strongly modified in case the QDs are put into 
the semiconductor microcavity (MC). In such structures 
electromagnetic wave is confined between two mirrors, 
what strongly modifies the light behavior for example 
stop-bands and standing electromagnetic waves can be 
observed^i^. 

Such systems are good objects for investigation of 
light-matter coupling. Depending on the interaction 
strength between QDs and MC two regimes are possible. 
First - weak-coupling regime shows slight modification of 
the photon dispersions. Photons and QDs excitons can 
be considered as different particles in this case. For large 
enough density of the quantum dots and small their in- 
homogeneous broadening one can expect also the called 
strong-coupling regime. It's characterized by strong in- 
teraction constant and it gives strong modification of the 
particles behavior. However in this paper we will con- 
sider the week coupling regime. 

Introducing to the system the acoustic fieldi"— one can 
observe significant effects, like acoustically driven ampli- 
fying of the lasingi. Acousto-optical interaction is based 
on a fact, that width of a band gap in semiconductors 
depends on lattice constant which value can be driven by 
acoustic pulse. And so the spectral possition of quantum 
dots could be modified by acoustic vibrations. Introduc- 
ing a strain pulse to the system of quantum dots coupled 
to a microcavity during lasing one can involve more dots 
in the process of photoluminescence and increase the las- 
ing power. 

In our work we apply the theoretical model which de- 
scribes an ensemble of quantum dots in optical micro- 
cavity. Here we consider quantum dots as a number of 



two-level systems, so the evolution of QDs population can 
be described by von Neumann equation. For the electro- 
magnetic field of the microcavity we can use a resonant 
mode approximationSr^iS, and get a system of two cou- 
pled equations which gives us a behavior of the system. 

In first part of the paper we describe in more details 
the model we propose and its possibilities and limitations. 
We show several results that could be obtained with this 
model. And after that we introduce the acousto-optic 
interaction to our system to describe theoretically exper- 
imental results^. In the final part of this work we dis- 
cuss a proposal of implementation of the surface acoustic 
waves to this system and describe effects of dynamical 
lasing patterns formation. 



II. THEORETICAL MODEL 

A. Quantum dots 

We use the density matrix formalism in a scalar 
approximationiiii^ where each dot can be described by 
its own 2x2 density matrix p. The dynamics of the sys- 
tem is given by von Neumannn equation: 



ihp = [Hp] 



(1) 



Diagonal elements of p give the probabilities to find quan- 
tum dot in ground or excited state and non-diagonal 
ones give the correlations between the ground and ex- 
cited states, that are responsible for the magnitude of 
the quantum dot polarization. In fact Hamiltonian and 
density matrix depend on in-plane dot position. However 
for the electromagnetic wave propagating in z-direction 
(that is the subject of our study) field is homogeneous 
through X and y and all the dots are in the same field. 
Thus we can consider the quantum dots ensemble to be 
homogeneous in the plane. Because we neglect all the in- 
teractions between the dots except the one from the elec- 
tromagnetic field the dynamics of each dot depends only 
on the local electromagnetic field and the resonant dot 
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transition energy. The distribution of the dots over the 
speetrum is accounted for by the spectral density func- 
tion of the dots nx{E'^), with the dot optical transition 
energy E"^. 

Next approximation of our model is that we consider 
only two types of interaction with the quantum dots: first 
is the interaction of quantum dot with resonant elec- 
tromagnetic field in the cavity that is explicitly put in 
Hamiltonian. The second is external pumping and relax- 
ation of quantum dots, that is introduced as an additional 
phenomenological terms in the equation for the density 
matrix of i-th type of quantum dots: 



(2) 



The Hamiltonian of i-ih type of quantum dot with reso- 
nant energy E^ can be written in the following form: 



H 







E' 



(3) 



d^ - it is the dipole moment matrix element between the 
ground and excited states of i-th quantum dot, Sq - Sam- 
plitude of the cavity electromagnetic field on quantum 
dots layer, E^ - energy of the first exciton state in the 
QD. The dipole matrix element of the dot can be ex- 
pressed via martix element of the interband currents of 

QD of j-type: c?* = Its value is defined only by 

quantum dot structure and wc take it as a parameter in 
our model, the ujq ~ Ei/h is a frequency of QD optical 
transition. 

The last term in Eq. ([2]) responsible for external pump- 
ing and relaxations in QDs reads as: 



22-^22)71 

Pll72 
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(4) 



If one neglect the interaction with the cavity mode the 
nonresonant pumping drives the system to some station- 
ary state with = P22 within the characteristic relax- 
ation time 1/71. As a first approximation we assume 
that is the same for all types of quantum dots. The 
non-diagonal elements relax to zero with characteristic 
time 1/72- All these values are taken as parameters of 
the model. In order to describe the spontaneous emission 
in QDs ensemble we add a white noise to the off-diagonal 
elements of Hamiltonian ^ . 

One more important parameter describing the quan- 
tum dots is their spectral density distribution Uxe- We 
take it as a gaussian distribution over energy: 

<s = noe (5) 

here Eqd is the center of QDs distribution and /S.Eqd 
is its width. 



B. Equations for electromagnetic wave 

Dynamics of electromagnetic field in the microcavity 
in the presence of the pumped QDs is a rather compli- 



cated problem. In this paper we study modification of 
the lasing thresholds by the deformation pulse and are 
interested in not very strong deviation of the pump in- 
tensity from the threshold one. Because of inevitable 
inhomogeneity of the QD spatial distribution there are 
always the preferential spots where the lasing starts first 
due to optimal gain conditions. These spots being of fi- 
nite size feeds all cavity modes but are most efficient for 
the modes slowly propagating along the cavity i.e. the 
modes with in-plane wave vector close to zero^. Thus in 
what follows we will treat only normal cavity mode as the 
one that starts first to emit the coherent light. It should 
be understood that this assumption becomes invalid for 
the pump strongly ahead of the threshold value. 

The dynamics of the electromagnetic field Eo of the 
normal cavity mode acting on the QDs can be found from 
the following oscillator-like equation: 



dt 



(6) 



Here huc = Re{fiWc) ~ ijc is the resonant photon en- 
ergy of the microcavity for zero inplane wavevector. The 
imaginary part of cavity resonance 7c results from the 
finite cavity lifetime. Second term ^jJx comes from the 
interaction with the resonant polarization Jx induced in 
QDs layer and coherent with the cavity electromagnetic 
field Sq. This polarization can be written as 



(7) 



For monochromatic wave with frequency uj we can solve 
Eq.®, find the p\i and substitute it into Eq.® to get: 



Jx — ^ 
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(8) 



In case of zero (or small) broadening of the QDs ensemble 



ujQ the same for all quantum dots) this equation 



gives the standard polariton splitting equation for homo- 
geneous exciton line. 



h{uj — UJc)h(uj — CJg + 172) 
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(9) 

For the opposite case of the broad QDs energy spectra 
that is the subject of our study we can calculate the mod- 
ification of the imaginary part of the cavity resonance due 
to pumped QDs as 
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(10) 

< then the amplitude of electromag- 



netic field will decay with time H/^^-^-^ and the system is 
in absorption regime. In the cases when this parameter 
is positive the amplitude of the field grows and it means 
that the system is in generation regime. So, analyzing the 
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"i^^^ we could analyze the generation-absorption transi- 
tions. On the fig. [T] we've plotted the dependence of 
7^^-^ on detuning between cavity resonance and center of 
quantum dots distribution for three values of p22- 
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FIG. 1: (Color online). Absorption coefficient dependence on 
detuning for different values of pumping. 

From the picture it's clear that lasing threshold de- 
pends on pumping power as well as on detuning. 



C. Acousto-optical interaction 

To describe the interaction of quantum dots ensemble 
with acoustical vibrations wc introduce a shift of the QD 
resonant energy proportional to the local deformation. 

-^^22(0 = E'' + Estrain{t) (11) 

Here E strain (t) is the energy shift of exciton levels caused 
by the strain. 

Now we have the equation describing the behavior 
and interactions between all three components: quantum 
dots, electromagnetic wave and acoustic vibrations. 



III. IMPLEMENTATION OF THE MODEL 

Let's first consider the time-evolution of our model sys- 
tem without any strain pulse We take P22 = P22 = 

0.7 

at starting point and we consider the detuning S = 
huJc — Eqz) = 12meV which corresponds to slightly above 
lasing threshold regime. Results of the calculations arc 
shown in Fig. [21 The black curve shows the time evo- 
lution of the photoluminescence intensity, and the red 
dotted curve shows P22 for quantum dots that are in the 
resonance with the microcavity. One can see that both 
electric field and the occupation P22 oscillate in time with 
the same frequency but shifted in a phase by |- . These os- 
cillations show periodical transitions of the system from 
lasing to absorption regime and vice versa. This effect 
comes from assumption of systems homogeneity and for 
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FIG. 2: (Color online). Black curve - relation of photolumi- 
nescence intensity to its maximum value versus time. Red 
dotted curve - time evolution of P22 

realistic non-homogeneous cases we should average in 
time the intensity coming from different lasing spots. 

Typical time dependencies of the electromagnetic field 
intensities for different detunings S = HuJc — Equ are 
shown on Fig. [3l The black curve shows the case when 
the cavity resonance coincides with the QDs distribution 
maximum. The red and blue ones are for the detun- 
ings slightly above and slightly below lasing threshold 
respectively. Electromagnetic radiation in the absorp- 
tion regime is defined by spontaneous radiation only and 
it is a narrow-band noise with small amplitude. 
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FIG. 3: (Color online). Intensities of electromagnetic fields 
calculated for different values of detunings. 

The application of the acoustic pulse change drastically 
the dynamics of the microcavity. First we took acoustic 
vibrations in a form of model harmonic oscillations. We 
put the detuning 6 ~ hwc ~ Eqd = l2meV, that cor- 
responds to slightly above lasing threshold case. The 
amplitude of vibrations we took equal to 1 meV and we 
calculated several curves for different frequencies of os- 
cillations. The results for 100, 50 and 20 GHz are shown 
on Fig. m From this picture one can see that even small 
vibrations can strongly amplify the signal from microcav- 
ity. This effect comes from the fact that in the presence of 
vibrations in the system, more quantum dots are involved 
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in lasing process. The characteristic time of microcavity 
light couphng to the resonant QDs is proportional to the 
dots density. The amplification of cavity field depends on 
vibration frequency. When these vibrations are too fre- 
quent, the quantum dots just don't have time to respond 
on external field. 
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FIG. 4: (Color online). Intensities of electromagnetic fields 
calculated for difi'erent frequencies of harmonic acoustic vi- 
brations. 

In order to compare the results of our approach with 
the experimental worki we took their strain pulses and 
put it into our model. The profile of the strain pulse 
is shown in Fig. [5] (a). Incident strain pulse comes in 
the moment of approximately 750 ps. Its starts from 
compressive deformation that increases the gap and the 
energy of QDs excitons. At the moment of 2100 ps the 
pulse reflected from the surface comes which starts with 
decompression part. In the experimental worki authors 
surprisingly observed that amplification of luminescence 
intensity from the reflected pulse is much stronger that 
from the incident one. To check whether it's possible to 
obtain this effect within our model, we put the real strain 
profile in our code and calculated the time evolution of 
the electromagnetic field amplitude. One can see it on 
fig. [5](b). And this effect was successfully reproduced. 

We propose the following explanation of the physics of 
this effect. One can notice that between two pulses there 
still remain small vibrations in the system. Their am- 
plitude is not larger than 1 meV but as we have shown 
in the beginning of this section even such small oscilla- 
tions could strongly modify electromagnetic field in cav- 
ity. From Fig. [5] (b) it is clear that between pulses the 
amplitude of the field is much larger than before the inci- 
dent pulse. So the initial conditions are "better" for am- 
plification before the second pulse than before the first 
one when the amplification factor of the refiected pulse 
could be less than of the incident. 

To check our assumptions we cut the vibrations be- 
tween two pulses and put the new strain profile in our 
model. One can see it in Fig. |6] (a). In the figure |6] (b) 
we've shown the ratio between the intensity in the pres- 
ence of acoustic pulses (lacit)) and the average intensity 
without them (/q) for two cases: with (black curve) and 
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FIG. 5: (a) Temporal profile of real strain, (b) Evolution of 
the electromagnetic fields amplitude in MC. 



without (red curve) the interpulse vibrations. 

It is clear that the real amplification from the reflected 
strain pulse is less than from the incident. 



IV. MODULATION OF THE LASING BY 
SURFACE STANDING ACOUSTIC WAVES 

Surface acoustic waves (SAW) are interesting objects 
to study light-matter coupling in semiconductor systems. 
Because of its long wavelengths they produce slowly 
changing in time (quasistationary for excitons in most 
planar systems) periodic potential. Recently SAW were 
used in studies of polariton condensates^^"— , photolu- 
minescence from quantum wells^i etc. We propose to 
implement SAW to obtain different lasing patterns in 
the ensemble of quantum dots coupled to microcavity. 
If we will consider two counterpropagating in xy-plane 
acoustic waves with the same amplitude we could obtain 
standing SAW (SSAW). Furthermore we can produce 
two-dimensional SSAW by the interference between two 
orthogonaly propagating one-dimensional SSAW. The 
Estrain in this case will be the function of planar co- 
ordinates and time and can be written in the form: 



Estram [x, y, t) = As AW siu {k^x) sin (kyy) sin (ujsAwt) 

(12) 
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Here Asaw - amplitude of exciton energy shifting by 
SAW (several mcV in order), kx,ky - wave- vectors of x— 
and y— SAW, ujsaw - frequency of SAW. 

It was shown in Sec. Ilia that for the same pump in- 
tensity, transition to the lasing regime strongly depends 
on the relative position between the cavity resonance and 
quantum dots distribution maximum. So for given pump 
power there is threshold value of detuning Sth that when 
1^1 l£ Sth there is a gain in the system. And it is clear 
that using strain provided by planar acoustic waves we 
can obtain some regions of cavity where conditions of 
lasing would be satisfied and in others would not. So we 
could observe patterns of bright spots. To derive the de- 
pendence of lasing intensity on planar coordinates we've 
considered a simple model based on a pair of kinetic equa- 
tion: 
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FIG. 6: (Color online) (a) Temporal profile of real strain with- 
out interpulse vibrations, (b) Amplification of the PL inten- 
sities with (black curve) and without (red curve) vibrations. 
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(14) 
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FIG. 7: (Color online) Evolution of photoluminescence in 
time. Snapshots are taken for the values of time 0.25, 0.52 
and 0.75 of the SAW period Tsaw 
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First equation describes number of photons (Nph) in 
system and second one - number of excited quantum dots 
{NxP22)- w - is a probability for the photon to be cap- 



tured by free quantum dot, Tc - lifetime of photon in 
MC when tqd is a non-radiative lifetime of an excited 
quantum dot. P is a term describing the pump and in 
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general case it should be proportional to the number of V. SUMMARY 

free quantum dots: P = pN^^l — ^22)- Here we consider 
only quantum dots which take participation in lasing, in 
other words which frequency of transition is equal to the 
MC frequency: 



QD 



Estraini,^-; ^)) 



(15) 

In the regions of absorption regime we can take Nph — 
0. In other regions we could derive Nph finding stationary 
solutions of Eqs.dlMll): 



ph 



2tqd 



PTQD + 1 

2wtqd 



(16) 



From this equation conditions for the lasing could be 
derived. First - pump should be strong enough so the 
prefactor before Nx to be positive: prgo > 1- Second, 
even with strong pumping it is necessary to have enough 
quantum dots to make Nph positive, so Tc(^^^^^)A^i: > 

PTQD+l 

2WTQ0 ' 

On Fig. [7] we plotted the solutions of Eqs. (jlSlllSp for 
different moments of time. We consider the pump cor- 
responding to the threshold detuning 6th = lOmeV, and 
the unstrained detuning of system we took to be slightly 
more. Asaw we took equal to 5meV. The video of time 
evolution of photolumincsccncc could be seen at^. 

Now let's consider weaker pumping to obtain threshold 
detuning Sth less than double strain amplitude. If we 
start with absolute value of detuning S slightly larger 
than 6th then the strain could be so high that \6\ will 
cross the threshold twice. In this case one can obtain 
holes in bright spots of photoluminescence as it is shown 
on fig. [H Video of this process could be seen at'^ . 



In our work we've developed theoretic model of light- 
matter interaction in the system of quantum dots coupled 
to microcavity in the presence of acoustic deformations. 
We've described the effect of acoustically driven ampli- 
fication of the lasing regime in the system. Moreover, 
we successfully reproduced experimental effects obtained 
iiii. At the final section of our work we proposed an im- 
plementation of surface acoustic waves to modify lasing 
patterns of the system. 



Acknowlegements 




FIG. 8: (Color online) Snapshop of PL for S ~ Sth = 2meV 
and Asaw = 5meV 



This work was supported by EU ITNs Spin-Optronics 
Grant No. 237252 and INDEX Grant No. 289968. 
Also we would like to thank C. Bruggemann, A. V. 
Schcrbakov, A. V. Akimov, D. R. Yakovlcv, M. Bayer, 
I. V. Ignatiev, V.D. Kulakovskii and G. Mapluech for 
useful discussions. 



7 



^ Zh. I. Alferov, Semiconductors 32 (1) (1998). 

^ A. V. Kavokin, J. J. Baumberg, G. Malpuech, F. P. Laussy, 

Microcavities (Oxford University Press, 2007). 
^ C. Weisbucii, M. Nishioka, A. Ishikawa, and Y. Arakawa, 

Phys. Rev. Lett. 69, 3314 (1992). 

* A. V. Akimov, A. V. Scherbakov, D. R. Yakovlev, C. T. 
Foxon, and M. Bayer, Phys. Rev. Lett. 97, 037401 (2006). 

^ A. V. Scherbakov, P. J. S. van Capel, A. V. Akimov, J. I. 

Dijkhuis, D. R. Yakovlev, T. Berstermann, and M. Bayer, 

Phys. Rev. Lett. 99, 057402 (2007). 
® T. Berstermann, A. V. Scherbakov, A. V. Akimov, D. R. 

Yakovlev, N. A. Gippius, B. A. Glavin, I. Sagnes, J. Bloch, 

and M. Bayer, Phys. Rev. B 80, 075301 (2009). 

C. Bruggemann et al.. Nature Photonics 6 (2012). 

* N. A. Gippius, T. Weiss, S. G. Tikhodeev, and H. Giessen, 
Optics Express 18, 7569 (2010). 

^ N. A. Gippius, S. G. Tikhodeev, V. D. Kulakovskii, D. N. 
Krizhanovskii, and A. I. Tartakovskii, Europhysics Letters 
67, 997 (2004). 

^° N. A. Gippius and S. G. Tikhodeev, Journal of Physics: 

Condensed Matter 16, S3653 (2004). 
" D. N. Krizhanovskii, G. Dasbach, A. A. Dremin, V. D. 

Kulakovskii, N. A. Gippius, M. Bayer, A. Forchel, Solid 

State Communications, 119, 435 (2001) 



V. D. Kulakovskii, D. N. Krizhanovskii, G. Dasbach, A. A. 
Dremin, N. A. Gippius, M. Bayer, and A. Forchel, phys. 
Stat. sol. (a) 190, 421 (2002) 

G. Ramon, U. Mizrahi, N. Akopian, S. Braitbart, D. Ger- 
shoni, T. L. Reinecke, B. D. Gerardot, and P. M. Petroff, 
Phys. Rev. B 73, 205330 (2006). 

E. A. Cerda-Mendez, D. N. Krizhanovskii, M. Wouters, R. 
Bradley, K. Biermann, K. Guda, R. Hey, P.V. Santos, D. 
Sarkar, and M. S. Skolnick, Phys. Rev. Lett. 105, 116402 
(2010). 

E A Cerda-Mendez, D N Krizhanovskii, K Biermann, R 
Hey, M S Skolnick and P V Santos, New J. Phys. 14, 
075011 (2012). 

Vishnevsky, D. V. and Solnyshkov, D. D. and Malpuech, 
G. and Gippius, N. A. and Shelykh, I. A., Phys. Rev. B 84, 
035312 (2011). 

Tetsuomi Sogawa, Haruki Sanada, Hideki Gotoh, Hiroshi 
Yamaguchi, and Paulo V. Santos, Phys. Rev. B 86, 035311 
(2012). 

See Supplemental Material at 



http://arxiv.org/src/1210.3199/anc for [sawl.avi] 



Supplemental 



Material 



at 



http://arxiv.org/src/1210.3199/anc for [saw2.avi 



